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calculate  SREMP. 

This  report  represents  one  area  of  the  PSR  research  effort  in 
SREMP.  This  document  was  prepared  as  one  volume  of  the  multivolume 
final  technical  report  for  DNA  under  contract  DNA  001 -85-C-0235.  The 
technical  monitor  was  MAJ  William  J.  Farmer. 
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SECTION  1 
INTRODUCTION 


Recently,  using  Green’s  function  techniques,  we  reduced  a  time- 
dependent  SREMP  problem  to  the  solution  of  a  set  of  Integral  equations 
for  tangential  values  of  electric  and  magnetic  fields,  E,  B,  on  the 
ground  surface  [Van  Alstine  and  Schlessinger ,  1986].*  Whea  substituted 
into  formal  solutions  for  E  and  B,  these  surface  values  give  a 
complete  solution  to  the  problem.  The  integral  equations  and  formal 
solutions  are  an  extension  to  three  dimensions  and  arbitrary  Compton 
currents,  of  earlier  results  using  Laplace  transform  techniques  for  E 
and  B  fields  depending  on  only  one  spatial  variable  (height)  generated 
by  Compton  currents  depending  only  on  time.  In  this  paper,  we  use  the 
new  three-dimensional  equations  to  generate  exact  solutions  for  spe¬ 
cial  cases  that  include  SREMP  fields  that  are  produced  in  three  dimen¬ 
sions  above  a  perfectly  conducting  ground,  SREMP  fields  that  are 
produced  in  three  dimensions  in  air  and  ground  when  air  and  ground 
conductivities  are  equal  but  time  varying,  and  SREMP  fields  in  one 
dimension  when  the  time-dependent  ground  conductivity  is  a  constant 
multiple  of  the  time-dependent  air  conductivity.  In  addition,  we  show 
how  the  new  three-dimensional  equations  contain  our  early  one¬ 
dimensional  results  in  two  independent  forms. 

First  we  simplify  our  formal  solutions  for  B  and  E  by  recasting 
them  in  a  new  form  that  makes  clear  their  connection  to  underlying 
vector  potentials.  We  then  solve  our  integral  equations  for  the  case 
of  an  infinitely  conductive  ground  and  use  the  solutions  to  obtain  B 
and  E  everywhere  for  that  case. 

Next,  we  solve  our  integral  equations  for  the  case  of  equal  but 
time-dependent  air  and  ground  conductivities  by  using  integral 
properties  of  the  Green’s  function.  We  show  that  the  exact  solution 


*Van  Alstine,  P.,  and  L.  Schlessinger,  Source  Region  Electro¬ 
magnetic  Effects  Phenomena,  Vol.  4,  New  Methods  for  Determination 
of  Three-Dimensional  SREMP  Environments,  Pacific-Sierra  Research 
Corporation,  Report  1588,  December  1986. 
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provided  by  our  method  is  that  given  by  the  full-space  vector- 
potential  (which  exists  only  in  this  case).  Finally,  we  examine  the 
(one-dimensional)  case  in  which  B,  E,  and  the  current  J  depend  only  on 
height  (or  depth)  and  time.  We  obtain  the  corresponding  formal  solu¬ 
tions  for  B  and  E  as  well  as  the  integral  equations  that  determine 
their  values  at  the  ground  surface.  We  show  that  when  J  is  further 
restricted  to  depend  on  time  only,  our  integral  equations  and  solu¬ 
tions  reduce  to  those  discovered  through  Laplace  transform  techniques 
in  our  early  work  [Schlessinger ,  198^].*  We  also  find  that  not  only 
do  the  one-dimensional  equations  and  solutions  emerge  as  the  remnants 
of  our  three-dimensional  equations  and  solutions  for  currents  and 
fields  that  only  depend  on  height  and  time,  but  because  of  the  struc¬ 
ture  of  Maxwell's  equations,  the  same  one-dimensional  equations  and 
solutions  govern  the  behavior  of  three-dimensional  fields  and  currents 
averaged  over  transverse  spatial  variables  (x  and  y). 

We  then  solve  our  new  one-dimensional  equations  for  height- 
dependent  J  for  the  case  in  which  the  (time-dependent)  ground  conduc¬ 
tivity  is  a  constant  multiple  of  the  (time-dependent)  air 
conductivity.  We  find  an  exact  integral  solution  that  is  general 
enough  to  include  previously  found  exact  solutions  for  the  cases  of 
equal  air  and  ground  (time-dependent)  conductivities,  time-dependent 
air  but  infinite  ground  conductivities,  and  unequal  but  constant  air 
and  ground  conductivities. 

Using  particular  realistic  forms  for  the  Compton  current,  we 
evaluate  our  solutions  for  the  fields  in  three  dimensions  above  a 
perfect  conductor  and  the  fields  in  one  dimension  when  time-dependent 
air  and  ground  conductivities  differ  by  a  constant  multiple.  These 
explicit  solutions  provide  insight  into  the  general  behavior  of  SREMP 
fields  as  well  as  analytic  test  solutions  for  comparison  with  numeri¬ 
cal  solutions  for  more  general  cases.  In  a  subsequent  report,  we  will 
present  the  numerical  results  for  SREMP  fields  obtained  from  these 
solutions  and  compare  them  with  those  obtained  by  other  methods. 

*Schlessinger ,  L. ,  Electromagnetic  Effects  Phenomena,  Vol .  1, 
Analytical  Solutions  for  SREMP  Environments,  Pacific-Sierra  Research 
Corporation,  Report  1437,  November  1984  (subsequently  published  by  the 
Defense  Nuclear  Agency,  Washington,  DC,  as  DNA  TR-84-397-V1 ) . 
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SECTION  2 

SIMPLIFICATION  OF  FORMAL  SOLUTION 


In  this  section,  we  use  vector  Identities  to  rewrite  the  integral 
expressions  for  B  and  E  ♦  J/o  of  Van  Alstine  and  Schlesslnger  [1986]* 
in  two  simpler  forms.  The  first  [Eqs.  (7)  and  (9)  below]  makes  clear 
the  role  played  by  image  currents  and  is  valid  everywhere.  The  second 
[Eqs.  (11)  and  (12)]  displays  B  and  E  ♦  J/a  as  curls  of  vector  quan¬ 
tities  (everywhere  but  on  the  ground  surface)  and  thus,  explicitly 
demonstrates  their  divergencelessness  (except  at  the  ground).  Both 
possess  simple  limits  for  infinitely  conductive  ground. 

In  Van  Alstine  and  Schlessinger  [1986], *  we  used  Green's  function 
techniques  to  derive  a  formal  solution  for  B  above  and  below  a  plane 
ground  surface  provided  only  that  the  (time-dependent)  air  and  ground 
conductivities  do  not  vary  in  space  and  that  the  displacement  current 
is  negligible  in  comparison  with  the  conduction  current  in  Maxwell's 
equations  with  Ohmic  conduction  current.  That  solution  was  given  by: 


B0, 


(1  -  2kk) 


dt’Q  •  (n  x  E) 


(D 


in  which  9g  is  8(+z)  and  n  »  -k  for  z  i  0  (air),  9g  is  8(-z)  and  n  - 
♦k  for  z  i  0  (ground),  Q  is  the  integro-differential  operator 


*Van  Alstine,  P.,  and  L.  Schlessinger,  Source  Region  Electro¬ 
magnetic  Effects  Phenomena,  Vol .  4,  New  Methods  for  Determination 
of  Three-Dimensional  SREMP  Environments,  Pacif ic-Sierra  Research 
Corporation,  Report  1588,  December  1^86. 

*Ibid. 
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where  erfc  is  the  complementary  error  function.  A  tilde  over  a  vector 
indicates  that  the  vector's  z  component  is  to  be  reversed  in  sign: 

v  *  (  V  v  ■  v  • 

The  identity 

(1  -  2kk)  •  (X  x  B)  -  A_  x  kB„  -  A  k  x  B_  -  kk  •  (Ar  x  B_)  ■  -  *  x  B  , 

T  Z  z  T  T  T 

(6) 


where  the  subscript  T  denotes  tangential  components,  allows  us  to 
rewrite  the  second  current  term  in  Eq.  (1)  as 


it 


-  (1 


-  2kk)  •  | V  X  Grj)  -  +  7  x y*^f  Gj  T 


Then,  Eq.  (1)  takes  the  simpler  form: 


B6B  -  -  7  x 


(Goj  -  GjT)  ♦  2jf  df 


a  •  (n  x  E)  , 


(7) 


which  makes  clear  the  role  of  the  second  current  term  as  the 
contribution  to  B  of  the  "image  current."  A  similar  rearrangement  of 
the  formal  solution  for  E  +  J/o, 


7’ 


x  J  - 


(1  -  2kk) 


x  J 


t 

f 


dt'  Q  •  (n  x  B) 


(8) 


results  in: 


x  J  -  Gj  7'  x  J)  -  2  /  dt'  Q  •  (n  x 


(9) 


t 


B) 


Since  the  operator  a  is  a  perfect  curl  everywhere  except  on  the  ground 
surface 


8  •  (n  x  E) 


•/ 


-7x7x|  dS'  n  x  EH 


(n  x  E)6(z)0(t  -  t' )  , 

(10) 


except  at  z  •  0  we  can  rewrite  Eqs.  (7)  and  (9)  in  terms  of  effective 
vector  potentials: 
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B6d  -  -  7  x 

D 


^  J (GoJ  "  Gi*^  +  2  7  XJ"  At'J ds'  n 


x  EH 


lz'-0 


(ID 


7*  x  J  *  Oj  V’  x  J) 


-  2  7  x 


/dt’f dS'  n  x  BHl 

J  Iz'-O 


(12) 


One  may  not  use  these  forms  to  match  solutions  on  the  ground 
surface  where  second  derivatives  of  H  are  singular.  Eqs.  (11)  and 
(12)  are  inaccurate  there;  they  should  actually  contain  the  6-function 
subtraction  included  in  Eq.  (10)  that  makes  their  formally  singular 
terms  equal  to  (the  nonsingular)  Q.  Therefore,  since  a  major  goal  of 
the  present  work  is  to  solve  the  integral  equations  for  B  and  E  at  the 
ground  surface  and  to  examine  behavior  of  the  fields  B  and  E  near  the 
surface,  we  shall  use  the  formal  solutions  Eqs.  (7)  and  (9)  that  are 
accurate  there  in  preference  to  Eqs.  (11)  and  (12). 


SECTION  3 

INFINITE  GROUND  CONDUCTIVITY  CASE 


In  the  case  of  an  infinitely  conductive  ground,  the  formal  solu¬ 
tion  for  B  leads  to  an  immediate  solution  of  the  Maxwell  equations. 
If  the  ground  conductivity  is  infinite,  the  electric  field  vanishes 
everywhere  below  the  surface.  Then  continuity  of  the  tangential 
electric  field  forces  the  tangential  electric  field  to  vanish  at  the 
surface  as  well : 


k  x  E  -  0  . 
z«0 


(13) 


Then,  the  electric  field-dependent  term  disappears  from  Eq.  (7)  giving 
B  in  the  air  directly  in  terms  of  the  Compton  current: 


-  V  x 


/ 


d4x' 

o' 


(GqJ  -  Gjft 


„  f  d4x’ 

J  ~ 


(10) 


where  Tg  «  1-jGd  ♦  kkGjj  is  the  "electric"  Green's  dyad,  while  Gp  and 
Gty  are  the  scalar  Dlrichlet  and  Neumann  Green's  functions, 
respectively.  Ampere's  law  (E  -  -J/o  ♦  V  x  B/yo)  then  leads  directly 
to  a  solution  for  E  in  the  air: 

GjX)  . 

(15) 


*As  in  Van  Alstine,  P.,  and  Schlessinger,  Source  Region  Electromagnetic 
Effects  Phenomena,  Vol.  4,  New  Methods  for  Determination  of  Three- 
Dimensional  SREMP  Environments,  Paclflc-Slerra  Research  Corporation, 
Report  1 ,  December  1986,  the  subscripts  >  and  <  denote  versions  of 
physical  quantities  evaluated  above  and  below  ground,  respectively. 
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But,  In  the  air. 


3^Gg  *  5(x  -  x* )6(t  -  t' ) 


(16) 


Equations  (14)  and  (19)  give  the  complete  solution  for  the  fields  B 
and  B  everywhere  above  an  infinitely  conductive  ground  in  terms  of  an 
arbitrary  Compton  current  distributed  in  the  air. 
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SECTION  4 

EQUAL  AIR  AND  GROUND  CONDUCTIVITY  CASE 


When  the  (time-dependent)  air  and  ground  conductivities  are 
equal.  Green's  theorem  applied  to  all  of  space  immediately  gives  the 
solution  for  8: 


>+< 


(20) 


and  hence  (through  Ampere's  law)  E: 


E 


1_ 

o 


I  d4x'  CQj  - 

>*< 


V_  f  d4x’ 

uo  J  o'  u0 
>+< 


V' 


J 


(21  ) 


(when  J  is  continuous  at  the  ground  surface).  Note  that  in  this  case, 
the  Helmholtz  theorem  implies  that  the  divergenceless  vector  field  B 
is  given  both  in  air  and  ground  as  the  curl  of  a  single  vector  poten¬ 
tial  which  may  be  taken  to  be: 


>+< 


Thus,  since  we  already  know  Its  solution,  this  case  serves  as  a  check 
on  the  validity  of  our  integral  solutions  and  integral  equations.  It 
provides  us  with  an  opportunity  to  use  integral  identities  (that  may 
be  important  in  the  solution  of  the  general  problem)  to  effectively 
solve  the  integral  equation  that  is  the  heart  of  our  method.  Once  we 
have  solved  the  integral  equation,  in  the  rest  of  this  section  we 
shall  demonstrate  that  our  method  produces  the  result  contained  in 
Eqs.  (20)  and  (21). 

When  we  set  the  transverse  components  of  B  in  air  and  ground  as 
given  by  Eq.  (7)  equal  to  each  other  at  the  ground  surface  in  order  to 
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satisfy  the  appropriate  boundary  condition  (assuming  the  absence  of 
surface  currents),  we  obtain  the  integral  equation: 


/  «’»  •  (k  »  «>.<.  ,-O.T 

—  m 

in  the  general  case.  But  when  o>  -  o<  -  c(t),  the  air  and  ground  a 
terms  become  equal  so  that: 


That  is,  when  we  explicitly  write  out  the  meanings  of  the  surface 
values  of  terms  on  the  right-hand  side  as  inherited  from  matching 
transverse  B  in  air  and  ground: 


t 


Mow,  we  are  supposed  to  use  the  boundary  information  contained  in 
Eq.  (25)  to  determine  the  electric  field  term  in  Eq.  (7),  the  formal 
solution  for  B.  There  are,  in  fact,  at  least  two  and  perhaps  three  ways 
to  do  this.  The  one  that  we  shall  use  here  (because  It  Is  tailor-made 
for  the  present  case)  uses  integral  identities  implied  by  Green’s 
theorem  to  turn  Eq.  (25)  directly  into  an  expression  for  the  electric 
field  term  in  Eq.  (7).  A  second  method  uses  a  special  order  of  in¬ 
tegral  identities,  differentiations,  and  limits  to  invert  the  scalar 
Green’s  functions  contained  in  a.  In  one  dimension  this  operation 
determines  the  surface  value  of  the  tangential  electric  field,  which 
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when  substituted  into  Eq.  (7)  directly  gives  the  solution  for  B. 
However,  in  three  dimensions,  it  leads  to  a  more  complicated  vector 
object  which  is  itself  sufficient  to  determine  the  electric  field  term 
of  Eq.  (7).  Or,  we  may  try  even  in  three  dimensions  to  solve  for  the 
surface  value  of  tangential  E  eventually  using  it  in  Eq.  (7)  to  give 
B. 

Here,  when  we  apply  the  integral  identities  as  detailed  in  Appen¬ 
dix  B  to  Eq.  (25),  we  find  that  they  move  the  Q  term  off  the  z  -  0 
surface  to  variable  z,  giving 


(26) 


which  is  just  what  we  need  to  determine  B.  Substitution  of  this  into 
the  transverse  part  of  the  solution  for  B  as  given  by  Eq.  (7)  implies 
that 


®>T 


V  x  J  (GqJ,  -  Gj^)  -  V  x  f  GjJ^  +  V  x  f  ~~  GqJ< 

>  ?  .  ?  <; 


so  that 


>T 


V  x 


/  V1  v 

>+< 


(27) 


(28) 


Having  found  Bp  over  all  space,  we  could  now  evaluate  it  on  the  ground 
surface  and  use  the  result  in  the  z  component  of  the  Integral  solution 
Eq.  (7)  to  find  Bz  over  all  space.  However,  the  fact  that  B  is  diver¬ 
genceless  everywhere  provides  us  with  a  method  to  use  Eq.  (28)  for  B-r 
to  determine  Bz  directly  without  intermediate  surface  limits.  That 
is , 
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0 


(29) 


V 


implies  that 


3  B. 
z  >z 


-  V, 


>T 


(30) 


so  that 


z 


But  Eq.  (28)  states  that  Bp  is  the  transverse  part  of  the  curl  of  a 
vector  field: 


B>t  -  (V  x  A)?  , 


(32) 


where 


>+< 


When  we  substitute  Eq.  (32)  into  Eq.  (31) 


00 

/• 


a>*  -  /  «  ’ T  •  "  X  *>T 


H 


V  (?T  X  Azk)  ♦  VT  ♦  (k3z  x  At)]  . 


(34) 


we  find  that  the  integrand  is  a  perfect  differential  in  z,  so  that 


m 


Z 


(35) 
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(when  J  is  continuous  at  the  ground  surface),  after  use  of  the  Green's 
function  equation  and  integrations  by  parts  in  space  and  time. 


SECTION  5 

REDUCTION  TO  ONE-DIMENSIONAL  SOLUTION 


In  earlier  work  [Schlessinger,  1984],*  Laplace  transform  tech¬ 
niques  were  used  to  develop  a  version  of  our  solution  method  applica¬ 
ble  to  "one-dimensional"  problems  in  which  the  (three-dimensional)  B 
and  E  fields  depend  only  on  the  single  spatial  variable  z.  Such 
solutions  are  physically  important  for  two  reasons.  First,  for  (time- 
dependent)  conductivities  and  Compton  currents  of  interest  there  may 
be  periods  of  time  for  which  the  the  fields  vary  in  the  horizontal 
spatial  coordinates  only  over  distances  much  larger  than  the  effective 
diffusion  length  (or  skin  depth): 

D(t,  f)  -  2(L)1/2  .  (41) 

For  such  periods,  the  one-dimensional  solutions  will  provide  good 
approximations  to  the  true  fields.  Second,  as  we  shall  show,  in  our 
case  the  Maxwell  equations  for  these  one-dimensional  B  and  E  fields 
are  identical  in  form  to  the  partial  differential  equations  obeyed  by 
B,  E,  and  J  averaged  over  the  transverse  spatial  variables.  Thus, 
even  when  Compton  currents  and  the  fields  they  produce  are  strongly 
varying  in  the  transverse  variables,  those  fields  have  an  average 
(over  the  transverse  variables)  behavior  that  is  predicted  by  the  one¬ 
dimensional  problem. 

When  the  Compton  current  J  and  the  B  and  E  fields  produced  by  it 
depend  only  on  z  and  t,  each  of  the  spatial  integrals  in  the  formal 
solution  for  B,  Eq.  (7),  becomes  an  integral  of  the  transverse  spatial 
dependence  of  the  Green's  function  over  all  x  and  y,  i.e., 


Schlessinger,  L.,  Electromagnetic  Effects  Phenomena,  Vol .  1, 
Analytical  Solutions  for  SREMP  Environments.  Pacific-Sierra  Research 
Corporation,  Report  1437,  November  1984  (subsequently  published  by  the 
Defense  Nuclear  Agency,  Washington,  DC,  as  DNA  TR-84-397-V1 ) . 
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B0_  -  -  V  x 

D 


x  j  \ 7- dz'^  y*ds'Go) j(2’’ tM  ■  (/ ds,Gi) 

y  dtj[(/ dS '  G0  j  n  x  E(z',  t') 

V  J(/dS.H)nxE(Z..t.)  I  . 


J(z* ,  t' ) 


x  E(z ' ,  t ' ) 


Since  the  Green's  function  factorizes, 


(-— f 
\  2/0.1 


(x  -  x') 


0(t  -  t' ) 


(z  -  z')‘ 


—  \e'~ 

4L 

ZMJ 

\  2  sol 

:XT  -  V 


e(t  -  t' )  . 

(U3) 


and  since 


4L 


the  transverse  spatial  integral  of  G  is  just 


'dS'Go  ■  G01  ■ 


where  G^  is  the  one-dimensional  diffusion  Green's  function  that 
depends  only  on  z  and  z'.  The  consequences  of  Eq.  (45)  are  that 


dS'GQf (z ' ,  t')  -  G01f(z',  t’)  , 
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(47) 


3  /dS'G.f  (z'  ,  t')  -  3  G„,f(z',  t’)  , 
z  #  0  z  01 


and 


WT  •  y*dS'[(n  x  E)H]z,_0  -  0  . 


(48) 


Thus,  In  Eq.  (7) 


/• 


fl  •  (n  x  E)  -  /  dS'[n  x  EGQ]  -  [n  x  EGQ1 ] 

z '  «0 


(49) 


z'  »0 


so  that  Eq.  (7)  becomes 


%  •  -  “i  *  jf^r  -  =t,J>  *  2  f  «'U  *  K01]  . 


(50) 


Since  the  right-hand  aide  of  Eq.  (50)  is  transverse,  an  immediate 
consequence  is  that  Bz  vanishes  everywhere.  Since  Jz  does  not  con¬ 
tribute  to  the  current  term  of  Eq.  (50),  B  can  be  rewritten  as 


B0r 


-  k3z  x 


ff¥ 


dz'GDlJT 


/ 


dt’Cn  x  EGQ  ] 


(51) 


z'  -0 


Similarly,  performing  the  transverse  spatial  integrals  in  the  integral 
equation  [Eq.  (23)]  that  results  from  matching  transverse  magnetic 
fields  in  air  and  ground  at  z  -  0  yields  the  integral  equation: 


t 

/ 


dt'  k  x  E|  Gn. (z-0,  z'»0)  -  -  k3,  x 

z'-0  u 


/dt’ 

>-< 


dz'GQ1 (z-0)JT 


(52) 
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Likewise,  if  we  perform  the  transverse  spatial  integrals  in  the  formal 
solution  for  E  +  J/o  in  Eq.  (9),  we  find 

(e  +  -  )  0Q  -  —  ~  x  /  /  ^7-  dz’Gn1k3  ,  x  J  -  2 
\  o  /  B  no  z  J  J  o'  D1  z' 

(53) 

which  immediately  implies  that  Ez  -  -Jz  /a  everywhere.  When  we  carry 
out  the  same  integrations  in  the  integral  equation  that  results  from 
matching  transverse  components  of  E  as  given  by  Eq.  (9)  at  the  ground 
surface,  we  find 


(54) 


The  formal  solutions  for  B  and  E  and  the  integral  equations  contained 
in  Eqs.  (51)  through  (5*0  are  just  those  that  we  would  have  obtained 
had  we  applied  one-dimensional  Green's  function  techniques  directly  to 
the  one-dimensional  problem.  They  emerge  here  whenever  the  transverse 
spatial  variations  of  J  and  hence  B  and  E  are  small  enough  over  dis¬ 
tances  on  the  order  of  the  diffusion  scale  that  the  transverse  depen¬ 
dence  of  the  Green's  function  can  be  integrated  away. 

However,  there  is  a  more  general  significance  to  these 
equations.  If  we  return  to  the  formal  solutions  for  B  and  E  in  the 
general  case  as  given  by  Eqs.  (7)  and  (9),  and  average  them  over  the 
transverse  spatial  variables  using  Eq.  (45)  and  its  consequences,  we 
find: 
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■. If * 


<B>0b  ■  "  x  /  /  J7-  dz'GD1<JT>  +  2  /  dt’[n  x  <E>G01]  *  (55) 


f 


and 


(<*>  *  ) 


0B  ’  uo 


■  k3z  X//  —  dZ'GD1k3z'  x  <J> 


t 

-  2 y  dt'Cn  x  <B 


<B>001] 


z'  -0 


(56) 


where  <  >  indicates  an  average  over  transverse  spatial  variables. 

These  equations  are  identical  in  form  to  Eqs.  (51)  and  (53),  with 
transverse  spatial  averages  of  B,  E,  and  J  replacing  their  one¬ 
dimensional  versions.  Thus,  the  one-dimensional  solutions  have  a  more 
general  physical  significance  as  transverse  spatial  averages  of  three- 
dimensional  solutions,  in  fact,  this  is  a  direct  consequence  of  the 
structure  of  Maxwell's  equations  for  our  SREMP  problem.  When  the 
displacement  current  can  be  neglected.  Maxwell's  equations  read: 


V  •  E  -  J  (57) 

7  x  E  ♦  B  -  0  (58) 

V  •  B  -  0  (59) 

V  x  B  »  p(J  ♦  oE)  .  (60) 


Their  transverse  spatial  averages  then  become: 

w  -  ¥  (61> 


k3  x  <E_>  +  <B>  -  0 

Z  1 


(62) 
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3  <B  >  -  0 
z  z 


(63) 


k3z  x  <Rr>  -  u(<J>  +  o<E>)  ,  (64) 

when  a  -  o(t)  is  not  a  function  of  spatial  variables.  Note  that 
Eqs.  (61)  through  (64)  are  identical  in  form  to  the  Maxwell  equations 
[Eqs .  (57)  through  ( 60 ) ]  themselves  with  averaged  B,  E,  J  replacing 
one  dimensional  versions  of  B,  E,  J  that  depend  only  on  z  and  t. 

Equations  (51)  through  (54)  are  more  general  than  those 
originally  obtained  through  the  use  of  Laplace  transform  techniques 
[Schlessinger ,  1984].*  That  work  assumed  that  the  Compton  current  J 
depends  only  on  the  time.  In  that  situation  one  can  perform  an  in¬ 
tegration  by  parts  in  the  current  term  of  Eq.  (51)  so  that  one  can  use 
the  fact  that  3/3zJ(t)  -  0.  Then  Eq.  (51)  becomes: 


/dt'  (E*  ?)]z,.0  • 


which  can  be  rewritten  using  Ampere's  law  as: 


B0b  -  ±  2 


/uo’  (3z'BTG0l)z,.0± 


Similarly,  Eq.  (53)  becomes: 


(E  * » K  -  -  ib/dt'(n  *  fcoi)z,_0 


Equation  (66)  for  B  is  equivalent  to  Eq.  ( 1 6 )  of  Schlessinger  [1984] 
(after  some  rearrangement  of  the  earlier  work).  When  evaluated  on  the 


Schlessinger,  L.,  Electromagnetic  Effects  Phenomena,  Vol .  1, 
Analytical  Solutions  for  SREMP  Environments,  Pacific-Sierra  Research 
Corporation,  Report  1437,  November  1984  (subsequently  published  by  the 
Defense  Nuclear  Agency,  Washington,  DC,  as  DNA  TR-84-397-V1 ) . 
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ground  surface  z  -  0,  Eqs.  (66)  and  (67)  reduce  directly  to  the  In¬ 
tegral  relations  Eqs.  (17)  through  (20)  of  Schlessinger  [198#].*  If 
one  carries  out  the  same  integ*atlon  by  parts  and  uses  3/3zJ(t)  -  0  in 
the  intep*al  equations  [Eqs.  (52)  and  (5#)],  one  finds: 


and 


(68) 


(69) 


respectively.  Equation  (69)  Is  identical  to  the  Integral  equation 
[Eq.  (25)]  of  Schlessinger  [I98#],t  while  continuity  of  transverse  E 
in  Ampere's  law: 


*r 


1  *  ^1  T>|  r| 

z  A  “< 

-  — 1 

lz-0-*-  u0>  | z-0  °>  |z-0  ‘I z-0- 

u0< 

z-0  °<  1 

z-o 

(70) 


converts  Eq.  (68)  into  the  integral  equation  [Eq.  (2#)]  of  Schles¬ 
singer  [198#].*  Consequently,  all  of  the  analytic  and  semianalytic 
results  of  that  previous  work  may  be  obtained  from  the  more  general 
one-dimensional  solutions  and  Integral  equations  of  the  present  work 
[Eqs.  (51)  through  (5#)]. 
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SECTION  6 

GENERALIZED  ONE-DIMENSIONAL  EXACT  SOLUTION 


In  order  to  obtain  a  solution  to  the  Maxwell  equations  [Eqs.  (57 
through  (60)]  In  air  and  ground  In  the  general  one-dimensional  case, 
one  must  first  solve  one  of  the  Integral  equations  [e.g.,  Eq.  (52)] 
for  a  surface  field  and  substitute  the  result  into  an  integral  solu¬ 
tion  [e.g.,  Eq.  (51)]  to  obtain  a  field  throughout  space.  In  general, 
this  requires  the  numerical  solution  of  the  integral  equation. 

However,  for  a  special  class  of  cases  that  include  the  one-dimensional 
infinite  ground  conductivity  case,  equal  air  and  ground  conductivity 
case,  and  the  constant  but  unequal  air  and  ground  conductivity  case, 
one  can  use  the  integral  equation  [Eq.  (52)]  to  obtain  an  analytic 
solution  for  B  (and  hence  for  E) .  In  the  general  case,  the  integral 
equation  for  the  surface  value  of  E  is: 


) 


dt'k  x  E 


Iz'-O 


G01(z-0,  z’-0)>+< 


k3  x 
z 


ft 


fhd2’G01(Z-0)JT 


(52) 


In  order  to  use  this  information  to  determine  B,  we  must  transform  it 
into  the  term 


2 


dt’[n  x  EG.  ] 

u  z'  -0 


(71 ) 


in  the  formal  one-dimensional  solution  for  B,  Eq.  (51).  We  can  do 
this  either  by  solving  Eq.  (52)  for  k  x  E | z— 0  by  inverting 
y*dt'Goi (z-0,  z’-0)>+<,  or  by  using  the  integral  identity  Eq.  (126)  of 
Appendix  B  (as  we  did  in  Sec.  H)  to  turn  the  left  side  of  Eq.  (52) 
into  Eq.  (71).  We  choose  the  second  method  here.  Note  that  the 
integral  identity  Eq.  (126)  of  Appendix  B  tells  us  how  to  perform 
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integrations  involving  products  of  two  Green's  functions  that  solve 
the  same  differential  equation.  But,  if  we  multiply  Eq.  (52)  by 
either  3zGoi  >  or  ^z^OK*  we  be  forced  to  deal  with  integrals 
Involving  both  G>  and  G<.  Our  integral  identities  will  apply  to  these 
only  if 


G01<(z-0,  z'-0)  -  oG01>(z-0,  z’-O) 


(72) 


where  a  is  a  constant  in  space  and  time.  This  occurs  when 


L 


< 


^/a2  , 


(73) 


which  implies  that 


2 

a 


o<(t)/o>(t) 


const.  , 


(74) 


for  all  t.  In  this  case,  Eq.  (52)  becomes: 


t 

’/ 


0  ♦  a)  /  dt'k  x  E  GQ1  (z-0,  z'-O) 
z'-O 


-  11m 


lim. 


-*0*  k3z  *//^:<lz'O0t>JT> 

> 

-°-  «,  «//  Sf  dZ'G01<JT<  ’ 


(75) 


whose  left  side  is  entirely  given  in  terms  of  G>. 

Before  we  can  use  our  integral  identities,  however,  since  they 
refer  to  Green's  functions  whose  common  z  argument  has  been  made  to 
vanish  from  the  same  side  of  the  plane  z  »  0,  we  must  rewrite  the  term 
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on  the  right  of  Eq.  (75)  that  depends  on  J<  in  terms  of  a  z-*0  limit 
from  above  (just  as  we  did  in  Appendix  B): 


litnz-*0-  k3z  x 


Itf 


dZ'G0KJT<  “  -  UV 


-  -  lim 


0.  kaz  k  i/f-VV- 

0.  kSs  k//^-o,< 


•>T<<-r> 


(76) 


Because  of  the  form  of  Eq.  (76),  we  are  faced  with  a  further  problem. 
Even  though  we  have  assumed  in  Eq.  (72)  that  G<  is  proportional  to  G> 
when  both  z  and  z’  have  been  set  to  zero,  this  is  no  longer  true  when 
as  in  Eq.  (76)  z'  is  variable.  In  fact, 


(77) 


However,  if  we  use  Eq.  (77)  in  Eq.  (76),  we  find  that  because  the  z 
variable  is  eventually  set  to  0,  G<  may  still  be  turned  into  G>  by 
scaling  z'  so  that  z'new  -  a  z'0i<j'  That  is. 
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(z  -  Z')2 


The  integral  equation  [Eq.  (75)]  has  now  become: 


(1  ♦  a)  /  dt'k  x  Ej  GQ1>(z-0,  z'-0) 

J  Iz'-O 


lim  (  -  k3  x 

z-*0*  1  z 


ff%  «'°01>JT>  -  V  11  // “'WkIv) 


(79) 


We  multiply  both  sides  by  -23Z.G>  |Zf ,q  •  integrate  over  ./'dt,/po,>,  and 
use  the  integral  identities  Eqs.  (129)  and  (128)  of  Appendix  B.  The 
net  result  is  that 


t 

(1  +  «)/  dt’[k  x  EG01>]z’-0 

—  • 

' 1,32 1  ,1I'GI'>JT>  * ;  k8z  *  fjir  <■*'<=!, >JT<('-r ) 


(80) 


Setting  z'-*-z'  in  the  second  term  on  the  right  side  and  dividing  both 
sides  by  (1  ♦  a)  then  converts  this  into: 


2U 


t 

J  dt'Ck  x  EG01>]z,_0  -  — —  k3z  x  JJ  —  dz'GI1;>JT> 

> 

*  jfrr^T  k3*  k  jj*f  as'WT<(r)  • 


an  expression  for  the  surface  E-field  term  in  the  formal  one¬ 
dimensional  solution  for  B,  Eq.  (51).  If  we  substitute  Eq.  (81  )  into 
Eq.  (51 ) ,  we  find: 

B>  ’  '  k3z  x  ff  ^r-  dz'Go,>Ji>  +  k3z  x  ff^az'° h>jt> 

'  (1  a)  k3z  X  ~  dZ'GI1>JT> 

-  -r~ - r  k3  x  /*/*—•  dz'G.  J  (~)  (82) 

a(1  +  a)  Z  J  J  01  >  T<  \  a  / 

in  which  the  two  image  terms  may  be  combined  to  give: 


B  •  -  k3  x 

>  z 


a(  1  *  a)  *  jj  ”  az'C0t>JT<  (j")  ' 


for  B  in  the  air.  The  analogous  procedure  applied  to  B  in  the  ground 
yields: 
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(84) 


By  using  Ampere's  law  (E  -  -J/o  +  V  x  B /\io) ,  we  compute  the  cor¬ 
responding  E  fields  in  air  and  ground: 


(85) 


(86) 


respectively.  The  one-dimensional  cases  in  which  the  ground  conduc¬ 
tivity  is  infinite,  the  (time-dependent)  air  and  ground  conductivities 
are  equal,  and  the  air  and  ground  conductivities  are  unequal  but 
constant  all  satisfy  a</< -  const,  and  hence  produce  B  and  E  fields 
given  by  Eqs.  (83)  through  (86).  These  solutions  generalize  to  time 
dependent  o's  (whose  ratio  is  constant)  and  J(z,t)  the  solutions  found 
in  Schlessinger  [1984]*  for  slowly- varying  conductivities  and  z- 
independent  J’s. 

*Schlessinger ,  L.,  Electromagnetic  Effects  Phenomena,  Vol .  1, 
Analytical  Solutions  for  SREMP  Environments,  Pacific-Sierra  Research 
Corporation,  Report  l4l7,  November  1984  (subsequently  published  by  the 
Defense  Nuclear  Agency,  Washington,  DC,  as  DNA  TR-84-397-V1 ) . 
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SECTION  7 

EXPLICIT  SOLUTIONS  FOR  SPECIALIZED  COMPTON  CURRENT  DISTRIBUTIONS 


In  order  to  investigate  important  features  of  the  SREMP  fields 
yielded  by  our  equations,  we  evaluate  them  for  two  Compton  current 
distributions  of  physical  interest,  one  in  the  one-dimensional  situa¬ 
tion  of  Sec.  6  and  one  in  the  space  above  an  infinitely  conductive 
ground  as  treated  in  Sec.  3.  First,  for  the  one-dimensional  case  in 
which  air  and  ground  conductivities  differ  by  a  constant  multiple,  we 
choose  a  current  pulse  that  exists  only  at  one  time,  that  is  constant 
in  z  in  the  air,  and  that  is  exponentially  attenuated  below  the  ground 
surface: 


J(z,t ) 


J  «(t) 


J<6U)ez/X 


z  >  0 


z  <  0 


(87) 


When  we  use  this  J  in  the  appropriate  solution  for  B  in  air,  Eq.  (83), 
we  can  immediately  perform  the  time  integrals  to  obtain: 


1 


k3z  |C0t><t)  * 


(rH) 


A 

J. 


5.1 

aO^-STT  k32  * / 01 


(88) 


in  which  L(t.t')  has  been  evaluated  by  the  6-function  at  t'  *  0  so 
that 


G«,  ( t )  =  G0(L(t,  t’))|t,.0 


1 


2/ifurn 


(z  -  z* )‘ 
4L(t) 


(89) 
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Since  the  air  current  was  assumed  to  be  constant  in  z,  the  two  terms 
in  Eq.  (88)  that  depend  on  the  air  current  have  integrands  that  are 
perfect  differentials  and  can  be  integrated  immediately  to  give: 


A 

»>  •  -  \  2  -VW1'-0’  l>  *  ft-rs  k  *  ^  Gn>(2'-°-  21 


aTr~ * -.t;:  *  2  K\  I  a2'Wt,e 


f 


\L 

ctA 


(90) 


But  Gj(z'-O)  -  Gq(z'-O)  so  that  the  integrated  terms  can  be  combined 
to  yield: 


2a  k  x  Gni  (z' -0,  t) 


>  1  +  a  a>  01 > 


—r\ — - — \ —  k  x  J  3  /  dz'Gri1s(t)eaX  . 

a(1  +  a)ax  <  z  J  01 > 


(91  ) 


Now,  the  integral  inside  the  z  derivative  in  the  ground  current  term 
in  Eq.  (91)  can  be  rewritten  as: 


(z  -  z'  )2  zj_ 

f  (  1  \  '  S  aA 

I  -  /  dz’ - I  e  e 

J  \  2/sr/ 


When  we  complete  the  square  in  the  argument  of  the  z' -dependent  ex¬ 
ponential  it  becomes: 


28 


29 


where 


z  S 

*  .2 

f(z)  m  e  e  erfc 


U’/a 


1/2 


2 —  ♦  -2- 

1/2  A 


Similar  steps  for  this  special  current  lead  to 


B<  -  9(C) 


( 1  ♦  o) 


z _ 


a  e 


k  x 


/\ 

$  -  ^ 

>  a2! 


A 

J<  1 

* k  *  r(zx 


K) 


\  f  (  clz  )  - 

1  -  -  \  f (-OZ) 

)  \ 

a  / 

for  B  in  the  ground; 


A 

-  ~  5(t )  -  ^  r~ — . 

°>  W0>  °  *  a)  Jo/^ 


(>■-§)  a- 


z 

J  |  ’  4L> 

1  J<T  |  e _ L  f(z) 


A  a. 


for  B  in  the  air;  and 


(97) 


(98) 
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-L  .  »  s(t)  .  »<i) 

0<  U0<  ( 


a  j _ a 

'  *  *’  ( °<^ 


A 

2—  S  -  -  d  e 

—  J>  2  tV 

g  J 


d. 

z _  > 

(*F-fe  [<’*«) «-» *  ('  -  =)rt— *])! 


for  E  in  the  ground. 

We  can  also  obtain  useful  solutions  for  the  field  in  the  space 
above  a  perfect  conductor  when  we  choose  a  simple  current  form.  For 
this  case,  we  choose  a  Compton  current  density  which  is  only  in  the 
radial  direction  and  is  a  shell  expanding  at  the  speed  of  light--that 


J(x,  t)  -  J(r ,  t)r  -  JqRq  ^  6  '2'  Ct"  exp  (-  r/X)  r  . 

r , 


We  also  assume  that  the  air  conductivity  is  given  by  the  simple  form: 


o(t)  -  oQ  exp  (-  t/tQ)  . 


(101  ) 


With  this  form  we  find  that 


L(t ,  t')  -  — j  [1  -  exp  -(t  -  t')/tQ] 


To  evaluate  the  magnetic  field  in  this  case  we  use  Eq.  (I1*). 
Integrating  by  parts  and  using  the  fact  that  curl  (J)  -  0  we  find 


B  -  -  2 


/^/dS'G0 


(k  x  J) 


Using  the  radial  nature  of  the  current  in  Eq.  (103)  we  can  perform  the 
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angular  integrals  to  obtain 


8  •  £  ♦/  ( 1 4*2?)  J‘»-  *■> ' 

—  CO  0 

(104) 

where  It  is  the  modified  Bessel  function  of  the  first  order.  Using 
the  form  for  J  given  in  Eq.  (100)  we  obtain 


uJqRq3  a  Ydx  exp 

B  ,  -  ♦/ 

2/irctpd 


-  1)  -  ctx/A)  exp 


/r2  +  c^t2x2  \ 
-  \  4d2il  I 


I rctxsine 


M 

(105) 


2 

2d  i 


where 


d 


2 


uo(t) 


l 


exp 


We  can  easily  obtain  the  solution  for  constant  air  conductivity  from 
Eq.  (105)  by  setting  to  -*  «  so  that  o(t)  -  oq  “  const.  We  find: 


“Vo3  , 

fu°o  1 

|  1  /2  a  /*  eXP  1 

Ctx  j 
\  1 

UOn  2  2.2  2 1 

0  r  +  c  t  x 

*  1 

2/ir  ct 

1  t  ) 

J  x(1  - 

.3/2  6Xp 
x) 

4t  (1  -  x) 

(106) 


where  the  argument  of  ^  is 


yo^rcxsine 
2(1  -  x) 
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Ampere's  Law  (e  ■  -  -  *  -  V  x  b)  then  insnediately  gives  us  the  E  field 
'  0  uo  ' 

corresponding  to  the  B  field  of  Eq.  (10M)  generated  by  a  radial  cur¬ 
rent  above  a  perfect  conductor: 


•  frpcoseig  -  e(psin8I0  -  rl ^ ) J  , 


(107) 


where  the  modified  Bessel  functions  IQ  and  I1  have  argument  (H£|iH§)t 
For  the  current  given  in  Eq.  (100)  and  conductivity  given  in 
Eq.  (101),  we  obtain  the  electric  field  corresponding  to  B  in 
Eq.  (105): 


where  and  I.  have  argument  /  rctx?in8  \  We  can  obtain  the  E-field 

\  2d  l  i 

solution  for  constant  conductivity  corresponding  to  Bc  by  setting  tQ-*» 
so  that  o(t)  -  Oq  in  Eq.  (108).  We  find 
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5(r  -  ct)e 


uJ0R0  ( uo0  \  C  dx  l  ctx  \  \ 

*  wTT  '  l  ^77^ e‘p  <-  t>  t 

•  [?COS9I0  -  J(9ln9I0  -  I,)]  . 


»”o  (r5  •  AVI 

Ht  (1  -  x) 


(109) 


I  uo0rcxsine  \ 

where  IQ  and  1^  have  arguments  ^  ~2~  -  "x)~  /  • 

In  a  subsequent  paper,  we  will  present  numerical  results  obtained 
from  all  of  these  solutions  along  with  a  comparison  of  those  numerical 
results  with  solutions  obtained  by  other  methods. 


APPENDIX  A 

CHECK  OF  INFINITE  CONDUCTIVITY  SOLUTION 


We  now  check  explicitly  that  our  procedure  has  actually  manufac¬ 
tured  a  solution  to  the  Maxwell’s  equations  and  boundary  conditions 
appropriate  to  the  infinite  ground  conductivity  case.  First,  the  fact 
that  Eq.  (14)  gives  B  as  a  perfect  curl  implies  that  B  is  divergence¬ 
less  everywhere  above  the  ground  surface: 


V  •  B> 


0  . 


(110) 


Then,  the  properties  of  the  Green's  functions  in  Eq.  (14)  show  im¬ 
mediately  that  B  solves  the  diffusion  equation  (with  time-dependent 
conductivity  and  appropriate  source)  everywhere  above  the  ground: 


(irr  “  3t)  B>  ■  * v  x  f  [5{x  -  • t,)j  ■ 6(x  ■  -  t')j] 

>  J> 

-  -  V  x  1 0 ( z )  —  -  e(-z)  —  »  -  V  x  —  for  z  >  0  . 

L  °>  °<  °> 


V  x  —  for  z  >  0  . 
°> 


Since  E  given  by  Eqs.  (18)  or  (19)  was  obtained  from  B  through 
Ampere's  law,  E  and  B  given  by  Eqs.  (19)  and  (14)  satisfy 


+  J_  V  x  B 

S  +  a>  “  W0> 


by  construction.  As  a  consequence,  E  and  B  are  also  related  by 

?  x  E>  -  7  x  Y  J  dV[G0J  -  Gjjb  -  V  x  J  [3tGQJ  -  a^J] 


-  B>  • 
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Finally,  we  may  calculate  the  divergence  of  Eq.  (19)  to  identify  the 
electric  charge  density  that  accompanies  our  fields.  We  find: 


(114) 


everywhere  above  the  ground.  [This  agrees  with  the  divergencelessness 
of  the  total  current  (Ohmic  plus  Compton)  implied  by  Ampere's  law  in 
our  approximation.] 

We  must  now  check  that  the  boundary  behavior  of  E  and  B  given  by 
Eqs.  (19)  and  (14)  is  consistent  with  our  problem.  First,  Eq.  (19) 
implies  that  tangential  E  is  given  by: 


(115) 


Then  the  fact  that  both  the  Dirichlet  Green's  function  and  its  tangen¬ 
tial  gradient  vanish  as  z-*0  implies  that 

El  -  0  .  (116) 

lz-0 


Ez,  on  the  other  hand,  is  given  by  a  superposition  of  the  Neumann 
Green’s  function  and  the  normal  derivative  of  the  Dinichlet  Green's 
function: 


Since  B  is  just  -7  x  E  so  that  B  vanishes  everywhere  in  the  ground, 
Eq.  (120)  implies  that  our  Bz  is  continuous  at  the  ground  surface  so 
that  B  is  divergenceless  everywhere.  Finally,  the  transverse  part  of 
B  [according  to  Eq.  (14)]  is  Just: 


37 


It 


BT  -  -  VT  X 


/d4x'  _  ,  f  d«x '  _  . 

o'  °NJzk  k3z  *  J  o'  °DJT  * 


0  22) 


so  that  at  the  surface 


I  z-0 


or 


'  */>  /Z-0 

“i.o  • 

>  Z-0 , T 


(123) 


(124) 


Thus,  B  is  entirely  transverse  with  double  its  free-space  value  at  the 
ground  surface.  Since  B  vanishes  everywhere  below  the  ground  surface, 
this  implies  a  surface  current  at  the  ground  surface  given  by: 

^  z-0 


(125) 


APPENDIX  B 

SOLUTION  OF  INTEGRAL  EQUATION  USING  INTEGRAL  IDENTITIES 


In  Appendix  A  of  Van  Alstine  and  Schlessinger  [1986],*  we  applied 
Green's  theorem  to  two  Green's  functions  to  obtain  the  integral  iden¬ 
tity: 


G.(x,  x") 


-  2 


/Sv/as'  [>z. 


G0(x.  x')G0(x 


\  x")l 

-lz'-0  + 


(126) 


The  subscript  z'-0+  indicates  that  the  surface  integral  is  to  be 
viewed  as  the  boundary  of  a  volume  integral  in  the  upper  half  space 
z'  2  0.  By  differentiating  this  identity  and  using  the  facts  that 
Go  *  Gq(x  -  x*)  while  Gj  -  Gj(x  -  xT* ) ,  we  find: 

VGjCx,  x")  -  -  V'GjCx,  x") 

’  2  /». h  [WX’  *'  )^”Go(x’  *  *">]  *  (127) 

L  J  z '  »0  + 

which  becomes 

?G!(X'  x">  *  "2  /^•/dS'  [3z'G0(x'  X,)V,G0(X*'  x">]  +  * 

(128) 

Because  on  the  surface  z"-0  Gx  -  G0  and  3zGj  -  9ZG0,  Eqs.  (126)  and 
(128)  become 


Van  Alstine,  P.,  and  L.  Schlessinger,  Source  Region  Electro¬ 
magnetic  Effects  Phenomena,  Vol .  4,  New  Methods  for  Determination 
of  Three-Dimensional  SREMP  Environments,  Pacific-Sierra  Research 
Corporation,  Report  1588,  December  1986. 
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respectively. 

How,  Eqs.  (129)  and  ( 1 30)  show  that  integration  over  time  and 
surface  of  32«Golz'-0+,  with  either  Gq  or  its  gradient  evaluated  with 
both  of  its  spatial  arguments  on  the  surface,  essentially  reproduces 
the  Green's  function  or  its  gradient  with  its  first  spatial  argument 
moved  off  the  surface.  Such  an  operation  would  take  the  a  term  in  the 
Integral  equation  Eq.  (25)  (whose  Green's  functions  are  fully 
evaluated  on  the  surface)  and  turn  it  directly  into  the  Q  term  in 
Eq.  (7)  (whose  first  spatial  argument  roams  over  the  upper  half¬ 
space).  If  we  write  Eq.  (25)  out  in  full,  we  see  that 


so  that  we  need  to  deal  with  two  further  complications  if  we  are  to 
use  our  integral  identities.  The  first  is  that  our  integral  iden¬ 
tities  involve  two  Green's  functions  both  of  which  are  boundary  values 
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from  above  while  the  second  term  on  the  right-hand  side  of  Eq.  (131) 
contains  a  Green's  function  evaluated  on  the  surface  from  below. 
Therefore,  we  rewrite  it: 


lira  _V  x 
z  —0 


/d4x'  _  .  ,  .  »  /*  dUx’  -  .  ,  , 

-p~  GgJ  “  li®  +v  x  I  — p- 

0  'J  z-0  /  0 


/ '[  °°J< 


(x’)  . 


Then,  Eq.  (131)  becomes: 


i  /  dt'  /  dS ' k  x  E(x')Gn(x,  x' ) I 

z-0  \J  Iz'-O 

-  Vt  •  J dS'k  x  E(x’)H(x,  x* )  | 

'  K-o*’  x/  ^  v>  -  u'z,0*?x/  ^  V<‘*> 


-  VTVT  •  I  dS ' k  x  E(x')H(x,  x* ) 


Secondly,  we  have  to  check  that  H  in  the  second  term  on  the  left-hand 
side  of  Eq.  (133)  obeys  its  own  version  of  Eq.  (129): 

'2/  ^  /dS’  U’V*'  **•  *’• 

M  J  Jzf»0  tz"«0 

m 

■  "  2  J  K,00(x'  -/  C0(xV  t’.  t"’> 

1  t”  Jz’-0  ,zn- 


a 

/dt  rt  * 

—PT  G  (x,  x”,  t,  t"’)|  -  H(x,  x",  t,  t")| 

u  lz"-0  I z"-0 


0  34) 
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In  short,  if  one  inspects  the  identities  [Eqs.  (126),  (128),  (129), 
(130),  (131*)].  one  sees  that  the  net  result  of  multiplying  each  term 
in  Eq.  (133)  &y  -2  3zGo(z"',  z)|z_q+  and  integrating  over ./'dt/po./'dS 
will  be  to  shift  the  first  z  argument  of  all  terms  from  z-0+  to  vari¬ 
able  z  and  to  change  VGq  into  and  ^Gq  into  1C  I  on  the  right-hand 
side.  That  is,  performing  this  operation  on  Eq.  (133)  yields: 


y\t. 


a  •  (k  X  E). 


-  - 
-  ?  xj  GJJ>  -7  xj  ~~  Gj J,,  C3f' ) 

>  >  Jt 


in  which  z  is  variable.  Setting  z' -*-z*  in  the  second  term  on  the 
right-hand  side  of  Eq.  (135),  we  see  that: 


2  /  dt'Q  •  (k  x  EL 


~  f  d4x'  _  f  d4x'  .  . 

-  ’v  —  ciJ>  -  v  - /  — 

>  <  J 


(136) 


Then,  the  transverse  part  of  the  identity  Eq.  (6)  tells  us  that 


(A  x  B)t  -  -  (A  x  B)t 


so  that  the  first  term  on  the  right-hand  side  becomes: 


Thus,  Eq.  (135)  may  be  rewritten  as 


h  °  e>t  -  •  (’  *j f  ^  cit>  *  *  v< ) 


(139) 
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ATTN:  R  SCHEPPE 

BOEING  MILITARY  AIRPLANE  CO 
ATTN:  C  SUTTER 

BOOZALLEN  &  HAMILTON,  INC 

ATTN:  TECHNICAL  UBRARY 

BOOZ  ALLEN  4  HAMILTON.  INC 
ATTN:  L  ALBRIGHT 

BOOZ  ALLEN  4  HAMILTON,  INC 
ATTN:  DDURGIN 

CALSPAN  CORP 

ATTN:  LIBRARY 

COMPUTER  SCIENCES  CORP 
ATTN:  ASCHIFF 

E  SYSTEMS.  INC 

ATTN:  MAIL  STOP  3 

E-SYSTEMS.  INC 

ATTN:  J  MOORE 

EG4G  WASH  ANALYTICAL  SVCS  CTR,  INC 
ATTN  A  BONHAM 
ATTN  J  GILES 

ELECTRO  MAGNETIC  APPLICATIONS.  INC 
ATTN:  D  MEREWETHER 

GENERAL  ELECTRIC  CO 

ATTN:  CHEWISON 

GENERAL  ELECTRIC  COMPANY 

ATTN:  J  PALCHEFSKY,  JR 

GENERAL  RESEARCH  CORP 
ATTN:  JWYSS 

GRUMMAN  AEROSPACE  CORP 

ATTN:  SMITH.  TECH  INFO  CNTR 

GTE  GOVERNMENT  SYSTEMS  CORPORATION 
ATTN:  TECH  LIBRARY 

HERCULES  OEF  ELECT  SYSTEMS.  INC 
ATTN:  R  LAZARCHUK 

HERCULES.  INC 

ATTN:  W  WOODRUFF 

HONEYWELL,  INC 

ATTN:  LIBRARY 

HONEYWELL,  INC 

ATTN:  SR4C  LIBRARY/T  CLARKIN 

I  IT  RESEARCH  INSTITUTE 
ATTN.  I  MINDEL 


INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  CLASSIFIED  LIBRARY 
ATTN:  TECH  INFO  SERVICES 

IRT  CORP 

ATTN:  B  WILLIAMS 
ATTN:  RW  STEWART 

ITT  TELECOMMUNICATIONS  CORP 
ATTN:  R  SCHWALLIE 

JAYCOR 

ATTN:  LIBRARY 

JOHNS  HOPKINS  UNIVERSITY 
ATTN  P  PARTRIDGE 

KAMAN  SCIENCES  CORP 
ATTN:  KLEE 

KAMAN  SCIENCES  CORP 

ATTN:  LIBRARY/B  KINSLOW 

KAMAN  SCIENCES  CORP 
ATTN  E CONRAD 

KAMAN  SCIENCES  CORPORATION 

ATTN:  TECHNICAL  UBRARY 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

KAMAN  TEMPO 

ATTN:  DASIAC 
ATTN:  R  RUTHERFORD 

LITTON  SYSTEMS.  INC 
ATTN:  EEUSTIS 

LITTON  SYSTEMS.  INC 

ATTN:  J  SKAGGS 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  PEREZ 

LOCKHEED  MISSILES  4  SPACE  CO.  INC 
ATTN:  TECH  INFO  CTR  D/COLL 

LTV  AEROSPACE  4  DEFENSE  COMPANY 
ATTN:  LIBRARY 

LUTECH.  INC 

ATTN:  TOMLIV 

MCDONNELL  DOUGLAS  CORP 

ATTN:  TECHNICAL  LIBRARY 

METATECH  CORPORATION 
ATTN:  WRADASKY 

MISSION  RESEARCH  CORP 
ATTN:  EMP  GROUP 

MISSION  RESEARCH  CORP 
ATTN:  JLUBELL 
ATTN:  JR  CURRY 
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MISSION  RESEARCH  CORP,  SAN  DIEGO 
ATTN:  V  VAN  LINT 

MITRE  CORPORATION 

ATTN:  M  FITZGERALD 

PACIFIC-SIERRA  RESEARCH  CORP 

ATTN:  H  BRODE,  CHAIRMAN  SAGE 
2  CYS  ATTN:  L  SCHLESSINGER 
2  CYS  ATTN:  P  VAN  ALSTINE 

PHOTOMETRICS,  INC 

ATTN:  I  L  KOFSKY 

PHYSICS  INTERNATIONAL  CO 

ATTN:  DOCUMENT  CONTROL 

R  4  D  ASSOCIATES 

ATTN:  DOCUMENT  CONTROL 

R  4  D  ASSOCIATES 

ATTN:  LIBRARY 

R  4  0  ASSOCIATES 

ATTN:  JP  CASTILLO 

RAYTHEON  CO 

ATTN:  H  FLESCHER 

RCA  CORPORATION 

ATTN:  G  BRUCKER 

RESEARCH  TRIANGLE  INSTITUTE 
ATTN:  M  SIMONS 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  B-l  DIV  TIC  (BAOB) 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  G  MORGAN 

S-CUBED 

ATTN:  A  WILSON 

SCIENCE  4  ENGRG  ASSOCIATES 
ATTN:  B  GAGE 
ATTN:  BZETLEN 
ATTN:  M  BORDEN 

SCIENCE  4  ENGRG  ASSOCIATES,  INC 
ATTN:  V  JONES 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  W  CHADSEY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  E  ODONNELL 
ATTN:  PJ  DOWLING 

SINGER  CO 

ATTN:  TECH  INFO  CENTER 

SRI  INTERNATIONAL 
ATTN:  DARNS 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  F  LEOPARD 


TEXAS  INSTRUMENTS.  INC 

ATTN:  TECHNICAL  LIBRARY 

TRW  INC 

ATTN:  LIBRARIAN 
TRW  INC 

ATTN:  A  R  CARLSON 
ATTN:  J  PENAR 

UNISYS  CORP-DEFENSE  SYSTEMS 

ATTN:  TECHNICAL  LIBRARY 
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